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We study bound states of abelian gauge theory in D = 1 + 1 dimensions using an equal-time, 
Poincare-covariant framework. The normalization of the linear confining potential is determined by 
a boundary condition in the solution of Gauss' law for the instantaneous A° field. As in the case 
of the Dirac equation, the norm of the relativistic fermion-antifermion (//) wave functions gives 
inclusive particle densities. However, while the Dirac spectrum is known to be continuous we find 
that regular // solutions exist only for discrete bound state masses. The // wave functions are 
consistent with the parton picture when the kinetic energy of the fermions is large compared to the 
binding potential. We verify that the electromagnetic form factors of the bound states are gauge 
invariant and calculate the parton distributions from the transition form factors in the Bjorken limit. 
For relativistic states we find a large sea contribution at low XBj ■ Since the potential is independent 
of the gauge coupling the bound states may serve as "Born terms" in a perturbative expansion, in 
analogy to the usual plane wave in and out states. 
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I. INTRODUCTION 

The hadron spectrum is simpler than one would expect. Deep inelastic scattering shows important contributions 
from sea quarks and gluons, yet qq mesons and qqq baryons are successfully classified in terms of only their 
valence quark degrees of freedom. Dynamical features such as masses and magnetic moments are consistent with the 
non-relativistic quark model, even though (light) quarks are known to be ultrarelativistic. 

The situation is made even more interesting by the fact that relativistic bound states are generally difficult to 
address in field theory. The perturbative expansion, which works well for non-relativistic atoms in QED, appears 
to be excluded for QCD hadrons, notwithstanding the apparent similarities of hadronic and atomic spectra. The 
velocities of atomic constituents scale with the fine-structure constant, w ~ a. Thus relativistic dynamics {v ~ 1) 
appears to imply strong coupling (a > 1). In the absence of a small expansion parameter analytic studies of relativistic 
bound states are generally model dependent. 

A small (perturbative) coupling can nevertheless be compatible with a relativistic bound state. A simple example 
is provided by an electron bound in a strong external field. The strength of the field oc Ze provides a new parameter: 
We may in principle have Ze of O (1) while e 1. To all orders in Ze the bound states are determined by diagrams 
where the electron interacts only with the external field. At this level the spectrum is exactly given by the Dirac 
equation. Self-interactions of the electron are O (e^) QED loop corrections, to be added perturbatively. 

The Dirac bound states do include effects of virtual e+e^ pairs which are created in the external field. This is 
made explicit by time-ordering the electron propagators. When the electron scatters into a negative energy state it 
propagates backward for some time, creating a "Z" diagram. The intermediate state then contains an additional 
e~^e~ pair. The Klein paradox |5] demonstrated early on that the Dirac wave function cannot be interpreted as a 
one-particle probability density. The paradox is elegantly resolved by considering pair production effects using the 
methods of quantum field theory |3]. 

The Dirac spectrum can be classified in terms of a single electron degree of freedom [J. This demonstrates that 
the virtual pairs present in relativistic bound states need not be explicitly seen in their spectra, in apparent analogy 
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to the situation for hadrons mentioned above. The Dirac states allow to examine this interesting phenomenon and to 
apply it to states that are bound by mutual interactions between the constituents rather than by an external potential 
[5]. Briefly, the argument is as follows. 

The energy component of an electron's four-momentum p = (p^ ,p) is unaffected by scattering in a time- independent 
(static) external potential A^{x). If p° > —me initially then ^ —\/p^ + throughout the scattering. Hence 
the ie prescription at the negative energy pole of the electron propagators in any Feynman diagram is irrelevant. In 
particular, we may use retarded electron propagators cx (p*^ -I- \/p^~-i-~mf -I- ie)~^ without affecting the positions of 
the bound-state poles in the Green function, i.e., the spectrum. The time-ordered dynamics, however, will look very 
different: the electron propagates forward only, eliminating the Z diagrams and the explicit e+e" pairs. Only single 
electron Fock states are present at any time, distributed according to the standard Dirac wave function. Note that 
the change in ie prescription is allowed only in the absence of loop corrections, i.e., only at the level of the Dirac 
states without O (e'^) perturbative corrections. 

It is useful to formulate the above observation in the language of Hamiltonian field theory. If H is the QED 
Hamiltonian, where A = (A°, A = 0) is given by the static external potential, the Dirac bound-state condition 



H\M) = M\M) 



(1.1) 



has no simple solution in terms of standard Fock states built on the perturbative vacuum |0) . This is because H \ 0) 7^ 0, 
i.e., the Hamiltonian creates pairs in the vacuum. Consequently the eigenstates |M) must have Fock components with 
an unlimited number of e^e~ pairs, as we already concluded above. In order to get retarded (rather than Feynman) 
electron propagators we need to change the boundary conditions, i.e., the vacuum state |0) — > |0)^ such that the 
electron field operator annihilates the "retarded vacuum" , 



tp{x) 



0. 



(1.2) 



This eliminates pair production from the vacuum, H \0) 
(apparent) single-electron state 



\M) 



0. The eigenvalue condition (|1.1| is then satisfied by the 

dxTp^t = 0,x) 



(p{x) 
X{x) 



|0)j 



provided the c-numbered wave function [(/?, x] ^ satisfies the Dirac equation 



[ - iV • 7 + e^^A^{x) + m] 



ip{x) 



M7° 



ip{x) 



(1.3) 



(1.4) 



Since the Dirac state ( 1.3 1 is built on the retarded vacuum it has both positive energy components &^ |0)^ and negative 
energy components d |0)^. It was observed previously [6^ that scattering amplitudes defined using retarded boundary 
conditions give inclusive cross sections. Analogously, the norm of the Dirac wave function may be viewed as an 
inclusive electron density. 



Generalizing the Dirac state (1.3) to a bound state formed by two fermions we get [S] 

\P) = j dxidx2^piit = 0,a;i)exp [iP ■ (a?i + a;2)/2] *(a?i - X2)'ip2it = 0,0:2) |0)^ , 



(1.5) 



where the wave function ^' now has 4x4 Dirac components. We assume that the fermions have distinct flavors, 
/ = 1, 2, and define the boundary conditions such that ipi{x) |0)jj = ipl{x) |0)^ — 0. The binding should be generated 
by the mutual interactions of the fermions in order for the state to be translationally invariant, i.e., to have a well- 
defined four-momentum P ~ (E,P). In analogy with the Dirac case we expect that relativistic binding at small 
(perturbative) gauge coupling can occur through an O (e") contribution to the A^ component of the gauge field. This 
is achieved by imposing a non- vanishing boundary condition on the solution of the equations of motion. 



|x| 



lim Ff^^F^^x) ^ ~2A\ 



(1.6) 



where A is a universal constant. Requiring lim|^|_j.oo(VA°)^ = A'' in Gauss' law adds an O (eA^) instantaneous 
interaction term to the Hamiltonian, H ^ H + Ha, which for neutral states is effectively [5j 



Ha = -- 



;A2 



f,f' 



dx dyip]:'iljf{t,x)\x ~ y\i\j),i\jj\t,y). 



(1.7) 
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and thus is leading compared to the O (e^) perturbative interactions. The condition of stationary time development 
at Cl(eA2), 



(li,,^Ht.)\P)=E\P) 



(1.8) 



is satisfied by the /1/2 state (1.5) provided its wave function satisfies the bound-state equation 

iV, • {7%, - \P ■ [7%, ^(x)\ + mi7°*(a;) - m2^{x)-{' ={E- V{x)\ *(a;), (1.9) 

with the linear potential 



(1.10) 



The wave functions ^ and the energy eigenvalues E that solve the bound-state equation (1.9) depend on the 3- 
momentum P . In this equal-time, Hamiltonian formalism boost invariance (as well as time translation invariance) is 
a dynamical symmetry. Because the field theory (here QED) is Poincare invariant and is solved at lowest order in 
the coupling e, with the Poincare invariant boundary condition (1.6), we expect the state (1.5) to be covariant. In [7] 
we verified this in Z? = 1 -t- 1 dimensions using the boost generator M"^. The boosted state satisfies the bound-state 
equation (1.9) with the appropriately shifted momentum, i.e.. 



(1.11) 



Also in £) = 3 -|- 1 dimensions the energy eigenvalues of (1.9) have the required dependence on the momentum [S] 



E 



(1.12) 



This indicates that the states are Poincare covariant also in four-dimensional space-time. 

In order to keep the discussion as simple as possible we use the example of an abelian gauge theory (QED). The 
relevant application of this method is to hadrons in QCD, where a linear potential is desired. As shown in [51 [7] the 
generalization to a non-abelian symmetry is straightforward. 

The present paper is a sequel to our study [7] of Poincare invariance. Here we examine other features of the wave 
functions in D = 1-1-1. They have several unusual properties, some of which are shared with the (previously known) 
solutions of the Dirac equation. The = 1 + 1 wave functions can be expressed in terms of confluent hypergeometric 
functions, which allows detailed analytic and numerical studies. We expect that the results in this paper shed some 
light also on the properties of the solutions of (1.9) in 13 = 3 -)- 1. 



In the next section we study the solutions of the Dirac equation ( 1.4 ) for a linear potential in D = 1 + 1 dimensions 



We show how, for nearly non-relativistic dynamics (m > Ze)^ the solutions agree with those of the corresponding 
Schrod inge r equation at small fermion separations but reflect pair pr odu ction at separations where y{x) > 2to. 



In Sec. Ill we find the analytic solutions of the // bound-state equation (1.9) in D = 1 -I- 1 dimensions. Similarly to 



the Dirac wave functions they are not square integrable since their norm tends to a constant at large |a;|. The // wave 
functions, however, are generally singular at V(x) = E zL P^. Solutions that are regular at these points^ exist only 
for discrete bound state masses [Tlj . This differs qualitatively from the Dirac equation, whose solutions are regular 
at all X, giving a continuous mass spectrum. In Sec. IV we show that the normalization of the wave function at a: = 



can be determined by requiring duality between the bound state {S distribution) and fermion-loop contributions to 
the imaginary parts of current propagators. A self-consistent normalization is obtained in all Lorentz frames and for 
all currents. For highly excited bound states the wave functions [at low \x\, hence small t^(a;)] turn into plane waves 
of positive energy fermions as expected in the parton model. In Se c. [V| we evaluate the electromagnetic form factors 
and show that they are gauge invariant (in any dimension). In Sec. VI we express Deep Inelastic Scattering in terms 
of transition form factors in the Bjorken limit. For relativistic states (small fermion masses) the D = 1 + 1 parton 
distributions grow large at small XBj, with XBjfixBj) oc log^(xBj). Our conclusions are given in Sec. VII 



When the fermion masses are unequal, mi ^ m2 in 
pi -S. 00. 



1.9\, the wave functions arc fully regular only in the Infinite Momentum Frame, 
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II. DIRAC EQUATION IN D = 1 + 1 

Some of the novel properties of the // states that we study in I? = 1 + 1 dimensions, notably the asymptotically 
constant norm of their wave functions, are shared by electrons bound in an external linear A'^ potential. Soon after 
Dirac first proposed his wave equation [9 it was realized that solutions of the Dirac equation generally cannot 

be normalized, in contrast to the case of the Schrodinger equation, where the requirement of a finite normalization 
integral leads to the quantization of the energy spectrum. Thus in [12] it was shown that the solutions of the Dirac 
equation in Z? = 1 + 1 have a constant norm as \x\ — >■ cxd for all potentials of the form A'^{x) = (or combinations 
thereof), where n ^ is any positive or negative integer. A similar result holds in Z? = 3 + 1 dimensions when A — 0, 
for central potentials A°(r) which are polynomials in r or in l/r, with the interesting exception of A'^(r) cx 1/r [12j . 
The Dirac energy spectrum is thus generally continuous, and the completeness relation involves a continuous set of 
eigenf unctions [13]. 

The constant norm of the Dirac wave function reflects pair production in a strong potential. The phenomenon is 
related to the Klein paradox [2], which requires a multiparticle framework for its resolution [3]. The pairs can be 
seen to arise from Z diagrams when the scattering of an electron in an external field is time-ordered. As already 
mentioned in the Introduction, in the case of static A^ potentials the same spectrum is obtained with retarded electron 
propagators ^ . Then electrons of both positive and negative energies propagate forward in time, and the bound states 
are described by single particle Dirac wave functions. Due to the retarded boundary conditions the norm of the Dirac 
wave function is, however, "inclusive" in nature, in analogy with the more familiar concept of inclusive scattering 
cross sections jS]. 



The analytic solution of the Dirac equation (1.4) in D = 1 + 1 for a linear potential was first given in [TT]. We include 
the derivation below for completeness and to introduce our notation. We also discuss some numerical properties of 
the solutions, which appear not to be widely known. 

A. General solution 

We use a standard 2-dimensional representation of the Dirac matrices in terms of Pauli matrices, 

7° = 0-3 , 7^=10-2, 7*^7^=75 = 0-1. (2.1) 

In QED2 the potential generated by a static source of charge Ze is V{x) — eA^{x) = \Ze^\x\. We use units where 
Ze^ = 1, hence 

V{x) = \\x\. (2.2) 

In the following all dimensionful quantities can be given their physical dimensions through multiplication by the 
appropriate power of Ze^ . 



The Dirac equation (1.4) implies 



- idxX = [M -V - m)(p, 

-id^ip = {M-V + m)x. (2.3) 

We choose the phases such that f{x) is real and xi^) is imaginary. This means that the solutions are characterized 
by two real parameters, e.g., (p{0) and —ixiO)- By adding/subtracting the equations at x and —x we may impose that 

Lp{x) 7] (p{-x), x{x) = -V x{-x), (2.4) 

where rj = ±1. This allows us to consider solutions in the region a; > only. Continuity at a; = requires x(0) = 



for 77 ~ +1, and (f{0) = for = —1. The equations (2.3) then ensure that 92,c/?(0) = for rj = 1 and dxxi^) = fo^' 
r] = —1. 



The two first-order equations (2.3) give rise to a second-order equation for (p{x) 



d'Mx) + + [(^^ - - "^'] = 0, (2.5) 

where e{x) = x/\x\ is the sign function. For x ^ 00 the term V'^(p{x) must be balanced by d'^(p{x), hence 

^p{x 00) ^ exp(±ia;^/4). (2.6) 
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From (2.3) follows that x(a;) has a similar asymptotic behavior. Since the norms |<p(a;)| and tend to constants for 

X ^ oo the Dirac wave functions are not normalizable |10H12j , unlike the solutions of the non-relativistic Schrodinger 
equation. As we shall see, the solutions have features which support the interpretation that their norm at large V{x) 
reflects virtual pair contributions. 

The coefRcient of dx^{x) in (2.5) is singular at M — V + m = 0. Assuming if{x) ^ {M — V + m)" a.s M — V + m ~> 



we find a = or 2. Hence the general solutions ip{x) and x(^) have no singularities at finite x. Since the wave 
functions are not square integrable there is no restriction on the eigenvalues M. In Sec. |IIB| we discuss how the 
discrete eigenvalues required by the Schrodinger equation emerge nevertheless in the non-relativistic domain (m 3> 1). 
In Sec. |II C| we show that any two solutions with different eigenvalues M are orthogonal. 



Since only the combination AI — V{x) appears in the Dirac equation (2.3) it is convenient to replace x by the 
variable^ 



(2.7) 



The Dirac equation then reads'^, in the region where M — V{x) > 0, 



We may combine if and x iiito the single complex function 

(j){a) = [ipicr) + x(a')]e*'^ and conversely (p{a) = Re[(j){a)e~'"^] , 
The second-order equation for (/'(cr) is then 



X(fT)= J Im [0(^7)6-'" 



2a dl4> Aia)d^4> - 2m^4> = 0, 



with the general solution 



(cr) ~ (ai + ibi) i^TO^, i 2ia) + (02 + ib2) \fo \Fy{\ — \ivr? , |, 2icr), 



where \F\ is the confluent hypergeometric function and the a^, hi are real constants. From (2.9) we flnd 



^}(a ^ 0) = ai -I- a2\fo + 0(0-), 
X[a 0) = ibi -I- ih2\f^ + 0[a). 



Matching the terms of O (Xj \fa) in (2.8) gives the relations 

hi = 2mai, 



02 = 2mhi. 



(2.8) 

(2.9) 
(2.10) 
(2.11) 

(2.12) 
(2.13) 



The general solution of the D — 1 + 1 Dirac equation (2.3) with the linear potential (2.2) of QED2 is thus given by 

im^ 1 



■>P{<t) = (p{a) -l-x(cr) 



{a + ib)iFi 



,2ia^ + {h + ia)2m e{M ~ V)y^ iFi(^- 



im, 



-.2ia 



exp(— icr), 



2'2 7 ' ' ' /v-L.y^ 2 '2' 

(2.14) 

where a and b are real parameters and a = (M — V)'^. and x a-re given by the real and imaginary parts of the r/is., 
respectively. The sign function e(M — V) ensures that the solution is valid also in the region where M — V{x) < 0, 
since -y/o" in (2.8) changes sign at = 0. This solution agrees with Eq. (14) of |11) . 



The solution (2.14) is valid for a; > 0. The wave functions for a; < are given by the symmetry relations (2.4) 



Continuity at a; = requires that the anti-symmetric (real or imaginary) part of the wave function vanishes at a: = 0, 
which also ensures that the x derivative of the symmetric part vanishes. The positions a = aa where cither the real 



^ We consider solutions only for a; > in the following. The parity condition (2.4i gives the solutions for x < 0. 
^ Here and in the following we use the shorthand notation 1^(0") = </3[(T(a;)], and similarly for x{^)- 
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or imaginary part of the rhs. in (2.141 vanishes determine the mass eigenvalues M = 
set cc = at ctq. The eigenvalues AI thus depend on the ratio a/b of the parameters. 
The asymptotic behavior of the wave function at large x, 



ctq, since continuity allows to 



lim Tpia) = 



V4 



a + ife 



3 — i7r/4 



r[i(i + zTO2)] r(i + iim2) V2 



(2a) 



72, 



^7r(m^-i)/4 



/2cr 



a + ib 



r(-iim2) r[i(i 



im2)] 



(2a)^ 



72, 



(2.15) 



oscillates with a ~ a; /4 in agreement with the general result (2.6). Hence the norm of the wave function tends to 
a constant at high x. Unlike in the non-relativistic limit (the Schrodinger equation), the parameters of the solution 
cannot be determined by a normalizability condition. 



B. The non-relativistic limit 



The fact that the eigenvalues M of the Dirac equation (2.3) depend on the ratio a/b of the parameters in (2.14) 
raises the question about the approach to the non-relativistic limit. In Z? = 1 + 1 dimensions for a fixed potential this 
limit is equivalent to taking m — ^ cx), scaling simultaneously coordinates and momenta appropriately. For a linear 
potential the scaling of the coordinate and the binding energy Ej, = M — m is 



Eh 



-1/3 



In the non-relativistic limit the Dirac equation reduces to the Schrodinger equation 

whose normalizable solutions are given by the Airy function, 

p{x) ^ N Ai[m^^^{x-2Eb)] {x>0), 



(2.16) 



(2.17) 



(2.18) 



with p{—x) — ±p{x). The solutions are differentiable at a; = only for discrete binding energies. How are the same 
eigenvalues selected in the Dirac equation at large m? 

II + E(, in 



It turns out that the two independent solutions in (2.14) become degenerate at large m. Setting M 



(2.7) and noting the scaling (2.16) we have 



ii2Eb-\x\)] +o(to-2/3^ 



(2.19) 



At large m we may use a stationary-phase approximation in the integral representation of the Hypergeometric functions 
in (|2.14|), giving 



V'(ct) = (1 + ^)(a + 6)0FTOl/3g7rmV2+«(m^-V4)^j [^73(1^1 _ 2Eb)] 



l + O 



(2.20) 



which agrees with the standard solution (2.18) of the Schrodinger equation. Consequently the non-relativistic limit 
of the general solution ( |2.14 ) does not depend on the ratio a/b. 

As seen from Fig. [l] the Dirac eigenvalues are insensitive to a/b already for to > 1, where they merge with the 
bound-state mass given by the Schrodinger equation. At large to there is a very narrow range of a/b where the (real or 
imaginary parts of the) two terms on the rhs. of (2.14) nearly cancel. Only in this range does the position of the zero, 
a = (To, depend on the precise value of a/b. For, e.g., to, = 2.5 this occurs for the ground state near a/b = —1.041, 
and the width of the interval in a/b that gives a continuum range of masses M is of O (lO~^). The approximate 
cancellation in (2.14) then makes the wave function grow very rapidly near a — ao, mimicking the exponential growth 
of the Schrodinger solutions for general M. For generic values of a/b and m > 1 the bound-state masses given by the 
Dirac equation agree with those of the normalizable solutions to the Schrodinger equation as indicated in Fig. [l] 

The issue of the approach to non-relativistic dynamics was also addressed in |13j . A measure of the relativistic 
effects was provided by the distance from the real axis of certain poles related to the Dirac eigenvalues, which was 
found to be ~ exp(— kto^), with ir < n < 2tt. In view of this it is interesting to note that the (typical) difference 



7 



Dirac (b=0) 
Dirac (a=0) 
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m 



FIG. 1: Ground state masses M as a function of the c onstit uent fermion mass m (in units of V Ze^). (a) The dark (light) 
blue solid curves correspond to the Dirac wave function (2.14 1 with 6 = (o = 0). The dashed red curve shows M = m + Et, 



with Eb the eigenvalue given by the Schrodinger equation (2.171 (which is reliable for m > 1 only), (b) Absolute values of the 
ground state mass differences plotted versus on a logarithmic scale. The mass difference between the Dirac solutions for 
6 = and a = (solid blue line) decreases exponentially with m^, much faster than the difference of the 6 = and Schrodinger 
masses (dashed red curve). 



between the Dirac eigenvalues obtained with different a/b decreases similarly with m (here n = tt), as shown by the 
solid (blue) line in Fig. [TJb). 

As expected, the Dirac wave function is similar to the Schrodinger one only for values of x such that V{x) ^ to, 
i.e., for weak binding. Thus the (upper component of) the Dirac wave function agrees with the Schrodinger wave 
function in the non-relativistic regime, as seen in Fig. [2] (where to = 2.5 and 6 = 0). Following its decrease to very 
small values the Dirac wave function begins to increase at a value of x where V{x) ~ M, and becomes O (1) again, 
initiating its asymptotic oscillations (2.6) when V{x) ~ 2M. This is because the wave function depends on x through 
the variable a = [M — V{x)]'^, which takes the same value at x = as at F = 2M. The start of the oscillations at a 
potential energy corresponding to pair production is indicative of the relation between the non- vanishing asymptotic 
norm and multiparticlc effects (the Z diagrams mentioned above). 



C. Orthogonality 



Two distinct solutions of the Dirac equation (2.3 1, = i^k Xk) and $f with eigenvalues Mk ^ Mg satisfy 



icri9^$fc + mcrg^fc = (Mfc - y)^k, 



(2.21) 



Multiplying the first equation by $J from the left and the second by $fe from the right and subtracting, we find 

- ^94$l(Tl$fc) = {Mk - M,)($|$fc). (2.22) 
In terms of the solution '\l){x) = '^{x) + x(x) in ( 2.14[ ) this is 

9, [Im(i^,* Vfe)] = (Mfe - Af,)Re (^,* V^) . (2.23) 



8 




FIG. 2: The upper component 'fi{x) of the Dirac wave function (continuous blue line) and the Schrodinger wave function p{x) of 
( |2.18[ ) (dashed red line). The fermion mass is m = 2.5 and the parameter 6 = in the analytic solution ( |2.14[ |. Both solutions 
are normalized to unity in the range < a; < 6. 



Recalling that ip{x) = rj ip{—x) is real and xi^) = ~^ x{~^) is imaginary both sides of (2.231 are odd functions of x 
if Vk = —fje- Since wave functions of opposite parity are trivially orthogonal it suffices to consider the case r]k — rji. 
Integrating both sides from x — to x — X and noting that the insertion at a; = of the Ihs. vanishes we find 



lm[riMx = X)] = (Mk - Me) / dxRe[rMx)] ■ 



X 



The leading term in the asymptotic limit of (2.151 is of the form 

il;{x ^ +oo) ^T]tP*{x-> ~oo) = Ccr™'/2g 



(2.24) 



(2.25) 



where a = — Mx + and the complex constant C depends on m as well as on a and b, which need not be 

the same for each bound-state level. For a sufficiently large X that the asymptotic form (2.25) of the wave functions 
applies we have 



where Cke — C^Cke^^^'^i^ -^h) , If (to simplify the discussion) wc choose X such that^ 

cos[(Mfc - Mt)X] = 



(2.24) gives, taking into account that the integrand on the rhs. is symmetric in x, 



dx Tie\_^*fiipk{x)\ 



2Ke{CM) sin[(Mfc - Mi)X] 



-X Mk - Mi 

We recognize the rhs. as a standard representation of the S distribution. Taking X — >■ oo we have 



dx^l'^k = 



dxRe[i:;Mx)] = Re(Cfc,)27r^(A4 - Mi). 



(2.26) 



(2.27) 



(2.28) 



(2.29) 



The freedom in choosing the parameter a/b in the general solution (2.14) allows bound-state solutions for a con- 
tinuous range of masses M. Since all these solutions are orthogonal a completeness sum must include them all, as 
shown in 1131 . The spectrum of the Dirac equation with a linear potential is thus continuous, similar to the case of 
non-interacting plane waves. 



* Relaxing this assumption leads to an extra term in Eq. l|2.28|l, which is singular at = Mi but docs not contribute to the final result. 
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III. SOLUTIONS OF THE TWO-FERMION BOUND-STATE EQUATION IN L» = 1 + 1 



In this Section we give the general solution for the wave functions 5* = e"^<i> that describe the bound states (1.5) 
of two fermions in D = 1 + 1 dimensions, 



\E,P) 



dxidx2 exp [^iP{xi + X2)]'i{}i{0,xi)e^'^^{xi — a;2)V'2(0, X2) |0)j 



(3.1) 



The 2-momentum of the bound state is denoted {E, P), and the Dirac matrices are defined as in (2.1 ). The bound-state 
equation (1.91 is then ^[7] 



id, {ai, $(a;)} - {d^ip) {ai,^x)} - \P [(Ti,^{x)] + mia^^x) - m2$(a;)a3 = [E ~ V{x)\ 



(3.2) 



The linear potential imposed by the boundary condition (1.6 1 is of the form (1.10 1. We take the coefficient of }^\x\ as 
our energy scale and thus have the same potential (2.2) as m the Dirac equation, now with x = xi — X2- 

We showed in [T' that E = P'^ + and worked out the P dependence of the wave function. It turned out to be 
convenient to introduce the variable 



cj={E-Vf~P'^ = - 2EV + 



(3.3) 



which in the rest frame {P = 0) coincides with the variable a defined in (2.7), which we used in the Dirac equation^. 
The expression for a suggests to define^ the "kinematical 2-momentum" 11^ = — eA'^, where P'^ — {E,P)^ is the 
bound-state momentum and eA^^ = (y, 0)^, 



n(a;) = {E-V{x),P) = (coshC, sinhC)v^. 



(3.4) 



Here the last equality holds when a = > 0. The variable a{x) — <j(—x) first decreases with increasing x > 0, from 
cr = Af^ down to (T = — at V{x) — E, and then increases with x, behaving asymptotically as cr ~ a;^/4. 

The common phase (p{x), which for convenience is extracted from the 2x2 wave function $ in the state (3.1), is 

2 



ip{x) = e{x){ml - ■ml){( - = heix)iml - ?77|) log 



M 

E + P 



E-V + P 
E-V - P 



(3.5) 



Here C is defined^ in (3.4), and ^ is the rapidity of the bound state, — {E + P)/M . 

Taking the complex conjugate of the bound-state equation and changing x — —x we see that ^{x) and $*(— x) 
satisfy the same equation. Consequently we may define solutions of definite parity 77 = ±1 by 



$''(a;) = $(.t) + 77 $*(-a;). 



$'?(_2;) ==77 $''*(a;). 



(3.6) 



In the following we first construct solutions for a; > and then complete them to the region a; < according to ( |3.6[ ) , 
requiring continuity at a; = 0. 

The general structure of a 2 x 2 wave function $(x) that satisfies (3.2) is [7] 



^>(a;) = <!)+ -{rniW (t)-'m2(j3W), 
a 

(j){x) = $o(2^) + $i(a;)cri. 



(3.7) 



where $0 and $1 are scalar functions of x. Inserting these expressions into the bound-state equation (3.2) and 
substituting the variable a of (3.3) for x using 



d,^-{E-V)da 



{x > 0), 



(3.8) 



^ The invariant mass M is now that of the // system. The present variable a is related to the variable s of [7] through cr = _A/^ — 2e{x)s{x). 
® Here the space component of 11 has the opposite sign compared to its definition in [?]• 

Notice that | |3.4[ | defines f only for cr > 0. The last expression in l |3.5| l can be taken as the definition of ip for cr < 0, and this is enough 

to make the wave function # well-defined. 
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we find that <I>o and $i satisfy 



1 - 



(mi — r7i2)^ 



1 - 



(mi + m2)^ 



$i(a). (3.9) 



The expUcit dependence on E and P has disappeared, which means that $0 and <i>i are the same functions of a in 
any frame. They are, however, P dependent when viewed as functions of x due to the relation (3.3) between a and 



X. The fuU 2x2 wave function <i> is expressed in terms $0 and $1 by (3.7). The relation between the wave function 



in a frame where the CM momentum is P to the rest frame (P = 0) wave function is given by 



<i>(a) 



-.iC/2$(P=0)(^), 



(3.10) 



with ( defined by (3.4 1 and cr > 



A. // solutions for mi — 1x12 = m 



We first consider the equal mass case, mi 



m2 



m. Then the phase (p = in (3.5) and the coupled equations 



(3.9) reduce to a second-order equation for $i((t), which has the form of a Coulomb wave equation, 

4m^' 



49^$i - 

The solut ion will oscillate asymptotically, <1>i(ct ^ 00) 
function (2.6). The general solution for $1 is 



1 



$1 = 0. 



(3.11) 



exp(±i(T/2), analogously to the behavior of the Dirac wave 



$i(a) 



icr/2 



[a iFi{l - im^ 2, icr) + bU{l- im'^,2, icr)] , 



where a and b are constants and U{a, /3, z) is the Confluent Hypergeometric function of the second kind. 
The behavior of the U function for small argument. 



U{a,2,z-> 0) 



1 1 



T(a) z 



-+0(logz). 



(3.12) 



(3.13) 



causes the 2x2 wave function $(x) in (3.7) to be singular at cr = if 6 7^ in (3.12): Then limo— >.o ^i(o') = 
— i6/r(l — im^) is non-vanishing, and the singular factor 1/ct is uncancelled in (3.7). Such a singularity at cr = 
prevents even a local normalizability of the wave function, and causes the orthogonality integrals ( 3.35| ) (Sec. IIIC 
below) to diverge. 

The // bound-state equation thus differs significantly from the Dirac equation (2.3), even though both wave 
functions are oscillatory at large |a;|. The general solution for the Dirac wave function is re gular for finite x and thus 
locally normalizable, whereas this is true for the // wave function only provided 6 = in (3.12). 

If we express the bound state mass as M — 2m + then in the limit of large fermion masses (m — ^ cxi) the binding 
energy Ei, and the coordinate x scale as in (2.16) (in the rest frame, P = 0). Substituting a ~ Am? + 2m{Ef, — \\x\) in 
the solution (3.12) and using a stationary-phase approximation in the integral representation of the Hypergeometric 
functions they turn into solutions of the non-relativistic Schrodinger equation, 

2X2/3 

m 



ae-'"/^ iFi{\^im^,2,ia) = — 

\m, 

^^/2c/(l-zm2,2,za) --(2m2)2/3_ 



Ai 



^my^'{\x\-2Et 



TT e 



ae 



{At [(im)i/3(|:r| - 2^6)] + I Bi \{\mf'Wx\ - 2Ef,)\ } 



(3.14) 



r(l — im^) 

up to O {mT^/^^ corrections. The result for the U function involves the non-normalizable Airy Bi function. 

In order to ensure local normalizability*, orthogonality of the lowest-order solutions as well as the correct behavior 
in the non-relativistic limit we set 6 = and thus consider 

f 1 



$i(cr) = N a e^"'''^ iFi{l ~ im^ ,2,ia) = 



A^sinh(7rm 



-ia/2 



'(1 







$0(0-) = -$i(cr) - 2iNe'"'/^ iFi{l - im^, 1, ia) = -$5(cr). 



(3.15) 



The requirement of local normalizability was previously used in |14| . 
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where we assumed the normahzation constant N to be real This makes $i((t) real for all a, as may be seen by a 
M — > 1 — u transformation of its integral representation. Correspondingly, $o(f) is purely imaginary according to 



(13J|). 

The asymptotic behavior for large |f7| is 



<i>l(cr —T' ±00) 

<i>o(CT -> ±00) 



2 N 



TT TO 

2 N 

TT TO 



TO^logdcrl) + argr(l + iTO2) [l + 0{a-^)] 



(3.16) 



- le" 



9{-a) 



COS 



- TO^ log(|(T|) + argr(l + im^) 



\l + 0(a- 



where 6'(— ct) = (= 1) for ct > (ct < 0). Due to the oscillatory behavior of the wave functions, the magnitude of N 
cannot be fixed by a normalization integral. In Sec. |IV| we show that the normalization of highly excited states may 
be determined using duality between the contributions of bound states and free fcrmions to current propagators. 



So far we neglected an e{x) factor in (3.81 and thus assumed^ a; > 0. When the solutions for x < are defined 



according to the parity constraint (3.6) the choice of phase indicated in (3.15) implies that 



$i(a; = 0) = (7? = -l), 



$o(a; = 0)=0 (7? = +l). 



(3.17) 



The latter conditions ensure the continuity at a; = of <&i(a;), ^o{x) and their derivatives. They also determine the 
discrete bound-state masses M b y the positions of the zeros of $o(c) and <i>i((T) for t] = ±1, respectively, through 
a{x = 0) = ]VP according to ^3. 



In Fig.[3]ja) we compare the symmetric (rj = 1) ground state mass as a function of the constituent fermion mass 
with the solutions of the Schrodinger equation ( 2.18[ ) (at the reduced mass to/2). As expected there is good agreement 
for TO > 1. The mass difference decreases strictly monotonously with to, as shown on a logarithmic scale in Fig.lsF' 




Miff) - M{NR) 




(b) 



FIG. 3: The ground state mass M as a function of t he fer mion mass m [in units of the coefficient of the linear potential (2.2|]. 
(a) The solid blue curve is for the // wa ve function (3.15 1. The dashed red curve shows M — 2m + Efj, with Ei, the eigenvalue 
given by the Schrodinger equation (2.171 (with reduced mass \m). (b) The difference of the ground state masses in (a), plotted 
versus m on a logarithmic scale. 



In Fig. [4] we compare the shape of the (P = 0) symmetric <I>i(a;) ground state wave function in (3.15) for to = 4 
with the corresponding Schrodinger wave function ( |2.18 1. The two wave functions are nearly equal at low x where 
V{x) <^ TO. The relativistic // wave function increases from small values in the intermediate x region to begin its 
asymptotic oscillations near V = 2M. ^i{x) is symmetric in the region < V{x) < 2M since it depends on x only 
via(T = [M-V{x)]^. 



^ The solutions actually take the form of ( |3.15| also for a: < when mi 
mapping a = cr{x) has a kink at x = 0. 



: m2, but the extension to negative x is still non-trivial as the 
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FIG. 4: The ground state // wave function $1(2;) in ( |3.15 \ for P = and m = 4 (solid blue line) compared to the non- 
relativistic Schrodinger wave function p{x) in (2.181 with the reduced mass m — 2 (dashed red line). The argument of "I>i is 
a = [M - V{x)]'^ with M = 8.4100 while the binding energy for p{x) is — 0.4043. Both wave functions are normalized to 
unity in the region < a; < 5. 



The spectrum of // states is shown in the form of "Regge trajectories" in Fig. [5|[a), where the square of the 
bound-state masses are shown as a function of the excitation number n. The trajectories of the symmetric and 
antisymmetric (77 = ±1) solutions are degenerate, and in the case of small constituent mass m — 0.1 almost linear. 
For m — 4 the trajectory initially has a smaller slope. 

Notice that the lowest antisymmetric (77 = —1) state has M = for any m, since $1(0" = 0) = 0. This state was 
not included in Fig.[5ja). Fig. [5jb) shows that the wave function ^a{x) of this state is essentially non- vanishing only 
in the relativistic region, V{x)^ 2m. In the non-relativistic limit the wave function thus tends to zero. 




FIG. 5: (a) Regge trajectories for // bound states with equal constituent masses, m = 0.1 (upper curve) and m — 4.0 
(lower curve). The excitation number n of the first five symmetric {rj = +1, n = 0, 2, . . . 8, blue lines) and antisymmetric 
{ri — —1, n — 1, 3, ... 9, dashed red lines) states are plotted versus — {2m)^. (b) The wave function $0(2;) of the ri — —1 
state with Af = [not included in (a)] for m = 0.1 and m = 4.0. The normalization is chosen arbitrarily such that <l>o = — 1 at 
the first minimum. 



The mass spectrum can be solved analytically for high excitations, since then a{x = 0) = is large at the 
positions x = 0, where the wave functions (3.151 must vanish according to (3.17), and their asymptotic expressions 
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(3.16) may be used. The zeros are thus approximately at 



rr? log(a-) + arg r(l + irr? 



7r(n+i) (r/-+l) 
^ I Tin [j] = —1) 

This may be combined into the asymptotic mass spectrum 

Ml = 7rn + log(7rn) - 2 arg T (l + iw?) + O (n^^) , 

where n is odd (even) for r\ = +1 (j] = —1). 

On the other hand, in the h mit o f small fermion mass, m 0, we have u~*™^(l — uY 



(3.18) 



(3.19) 



in the integral representation (|3.15| for $i. Hence we find to O (m^) 



1 + im^ log[(l — u)/u] 



$i(ct) = iV(l - 



2 sin 



irn?a e 



'log 



/o ^ " 

For <T > the integral can be expressed in terms of sine and cosine integral functions as 



(3.20) 



-ia/2 



1 



u 



2 cos (I) [Ci(a) - log(a) - 7^] + 2sin (|)si((7), (3.21) 



where = 0.577216 is Euler's constant, and 

Si(z) 



2 



°° , sin(w) 
au , 



Ci(z) 



du 



cos(u) 



Imposing the continuity conditions (3.17) at x = we find the spectrum 



Ml^Tm + 2m2 [log(7r7i) - Ci(7rn) + 7^] + O (m"*) 



71 = 0,1,2, 



(3.22) 



(3.23) 



where n is odd (even) for 77 = +1 (77 = —1). The case rt = should be understood by taking the limit 71 — >■ 0, which 
gives Mq = 0, the solution shown in Fig. ^b), which is exact for any m. 



For TO exactly equal to zero the full wave function (3.7) reduces to <&(cr) — —2iN exp(z gicr/2 ), which is regular 



at all a. Hence there is no constraint on the spectrum when to = 0. On the other hand, (3.23) gives M„ = nn in 



the TO — !■ limit. The discrete spectrum obtained for regular solutions when m ^ thus differs, even in the to — >■ 
limit, from the continuous spectrum found with to = 0. Furthermore, the original bound state equation (3.2) implied 
parity doubling when TOi = TO2 = 0: The parity transformed wave function 7°<i>(— a;)7'' is a solution (with P —P) 
having the same eigenvalue E as $(a;). To the contrary, the to — >■ states have parity rj = (—1)"+^ and are not parity 
degenerate. 



B. // solutions for mi 7^ m2 



The general solution of (3.9 1 when 



Ato^ = — ni^ 7^ 



(3.24) 



IS 



$i(a) + $o(ct) = e"''^{\a\-i^"'\a+ib)miiFi[iml,l^il^TT? ,~ia) ~ |a|5'^™'(a-i6)TO2 iFi(wn?, l+^A7n^ -ia)}, 

(3.25) 

$i(cr) - $o(ct) = e''"''^{\a\^^'^\a-ih)mi iFi (-im^, I+^Ato^, icr) - |cr|-5^"'(a+i&)TO2 iFi(-iTO?, l-iATO^icr)}, 

where a and b are complex constants^°. The parametrization of the constants was chosen such that for real a and 6, 
<I>i is real and $0 is imaginary as in the case of equal masses above. We may assume that x > 0, with the solutions 
of definite parity defined as in (3.6). 



^'^ The absolute values in |(t|^*^™' /'^ specify the branch choices for positive and negative u. This choice is natural in view of our definition 
of the phase ip in K^, which behaves as e^f ~ |o-|=FiAm2/2 fo^ E - ^ ±P. 
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The 2x2 wave function $ is generally singular at cr = according to (3.7), which may be expressed as 



(3.26) 



Noting that iFi(a,/3,0) = 1 we find that for a =[{E -V) + P][{E - V) - P] ^ 0, the most singular terms of the 
full bound state wave function in (3.1) are 



(3.27) 



No choice of the parameters a, 6 can eliminate the 1/ct singularity at both E — V — P and i? — F = —P. The 
phase |f7|T*^™ , however, ensures the integrability of the wave function when multiplied by any regular function. For 
example, the orthogonality relations to be discussed below are well-defined. 

It is interesting to note that in the infinite momentum frame, P — > +oo (IMF), the full wave function is regular if 
a + lb = Q. This choice removes the l/ cr sing ularity in (3.27) at E ~ V — P, while in the IMF E — V —P at finite 
X. As P ^ oo both square brackets in (3.26) approach E{a3 + i(J2) = £'7+ and the full 2x2 wave function becomes, 
for b ~ ia in (3.25), 



e'^$P^oo(CT) = -2m AP^"' , £7+ e-'"/^ iPi(l - im^, 2 + iAm^,ia), 

1 + iAm^ 



(3.28) 



which is indeed regular at cr = 0, and also square integrable in x. While b — ia thus appears to be the most physical 
choice of parameters, we anyhow continue the discussion assuming generic a and b. 

By using known identities for the iPi functions it is straightforward to check that in the equal-mass limit mi, TO2 — ^ 
m the wave function $1(0-) in (3.25) reduces to 



$i(mi 



m2 



m) ibm [e'"''^ iFi{im^,l,-i(T)-e-'"''^ iFi{-im'^Xi<y)\ = -bma e''"^^ iFi{l-im^ ,2,ia), (3.29) 



which agrees with our previous expression (3.151 when N — —bm. The mi — >■ limit is also simple since iPi(a 
0,/?,z) = 1. Thus 



$j(mi = 0) = im2 [ - (a - ib) |a|-5"2eW2 ^_ _^ |^| 



(j = o,i). 



(3.30) 



The definition (3.6) of $(x < 0) requires continuity of ^oi^ — 0) and $1(0; — 0) for the bound-state equation (3.9) 
to be satisfied at all x, 



<^>]{x = 0) = 77$f(2; = 0) 
d^<i>]ix = 0) ^ -7?9^$f (a; = 0) 



U = 0, 1). 



(3.31) 
(3.32) 



The condition (3.31) requires that both $o(0) and $i(0) are real (imaginary) for rj — +1 {rj — —I). In general, this 



can be satisfied by adjusting the overall phase in (3.25) provided the phase difference of $o(0) and $i(0) is or n. 
This constraint determines the mass spectrum for fixed a/b. The continuity of the derivatives (3.32) follows from 
(3.31) and the fact that the wave functions satisfy the bound-state equation (3.9 1 for a; > 0. 



As mi — m2 the phase difference between <I>o and 'I'l approaches ±7r/2 as seen from Eqs. (3.9 ) and (3.29 1. Therefore, 



in the equal-mass case these wave functions cannot have the same phase, instead one of them has to vanish as we 



found in (3.17). As pointed out above, the same happens if a and b are both real (or, more generally, if they have 
the same phase, i.e., a/b is real). For mi « m2 the wave functions have a relative phase of nearly ±7r/2 everywhere 
except close to a zero of one of them. Thus the spectrum depends smoothly on mi — m2. 

In Fig.|6]we illustrate some properties of the solutions in terms of the density |<&oP + for the choice b — ia. The 
fcrmion masses mi = 1.0, m2 = 1.5 are fairly high compared to V'(x) = 0.5, ensuring that the multipair contributions 
to the wave functions are well separated in the ground state rest frame [solid red curve in (a)]. The wave function of 
the excited state [dashed blue curve in (a)] extends to larger fermion separations x before decreasing to small values, 
but its multipair contributions shift correspondingly in x, approximately preserving the extent of the gap. The same 
densities are shown in (b) for the case of non- vanishing center-of-mass momentum, P — 5. The wave functions Lorentz 
contract at low x, whereas the length of the gap to the multiparticle contributions grows. 
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P = 5 



M=5.11 



M=3.15 




10 



15 
(b) 



20 



25 X 30 



FIG. 6: (a) The density |$oP + |$i|^ as a function of the distance x between the constituents for the ground state (M — 3.15, 
solid red line) and for an excited state (M = 5.11, dashed blue line). The constituent masses are mi — 1.0 and m2 = 1.5. 
(b) The densities in (a) plotted in the case of non-vanishing center-of-mass momentum, P = 5.0. The densities are symmetric 
under x — ^ —x and normalized to unity at a; = 0. 



C. Orthogonality of the // states 



If we include the phase (p in the definition of the 2x2 wave function in the state (3.1 1, 

^{x) EE e*'^(^)$(a:), 

the inner product of two states k, £ reduces to 



{Ei, P, Pk) = 2tt5{Pu - Pi) J dxTv (x)^'fc(x)] 



since only the anticommutators of the fields contribute due to (1.2). Thus orthogonality of the states requires 



J dx Tr [^l (x)] =0 iPk=Pi^ P, Ek + El 



(3.33) 



(3.34) 



(3.35) 



The bound-state equations for (a;) and '^\{x) are 

i9^{CTl,*fe}-iP[(Tl,*fc]+TOl(73 4'fe-m2*fcCr3 = [Ek''V]-^k: 



id.,[a^,^l^ + \P 



0-1,^' 



(3.36) 



Multiplying the first equation by from the left, the second by from the right and then taking the trace of their 
difference gives 



(3.37) 



Integrating both sides over all x the Ihs. gets a contribution only from x — ±oo [only the $o and $i components of 
the 2x2 wave function $ in (3.7) contribute on the Ihs.]. To leading order in the x — -l-oo limit we have from (3.5) 
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and (3.251 



tp{x — > oo) 



The result in the 



Am^ log[M/(£; + P)], 

1 [1 + im\) 1 (1 + iTTij) 



C2 — m2{a — ib)- 



r (1 + iArr??) 
r(l - iml) 



e™?/2 ~mi(a + i6) 



r (1 - iAm^) 
r(l - im\) 



2/2 



(3.38) 
(3.39) 

(3.40) 



-00 h mit is given by the complex conjugate of the above, since "^{—x) — 77^'*(2:) according to 



(3.6). The product ^g^k in (3.37) oscillates asymptotically, and its integral may be defined analogously to that of 
plane waves, e.3., by adding a factor e"*^'^' with infinitesimal e > 0. Then the integral of the Ihs. of (3.37) vanishes, 
and the orthogonality (13. 35^ is ensured. 



IV. DUALITY 



A. Wave function normalization 



The bound-state equation does not determine the overall normalization (or phase) of the wave functions. Due 
to contributions from an infinite number of particle pairs the integral of the norm of the Dirac-type wave functions 
diverges. The relative normalization of high-mass states is needed to determine the parton distributions of the bound 
states in Sec. |VI| Here we shall use an approximate duality relation to determine their normalization. For simplicity 
we limit ourselves to the equal-mass case in the following, nii — m2 — m. 



k . 




k-P - 

FIG. 7: Duality between resonance and fermion loop contributions to the imaginary part of a current propagator. The relation 
should hold in a semi-local sense, and become more accurate at high excitations. 

In our present approximation the bound-state spectrum is a sequence of zero-width resonances. The (properly 
averaged) bound-state contribution to the imaginary part of a current propagator is expected to equal the contribution 
of a free fermion loop, as illustrated in Fig. [7] This allows to express the wave function at the origin in terms of 
a calculable perturbative loop contribution. We get consistent normalizations in all frames for scalar, pseudoscalar, 
vector and pseudovector currents. Furthermore, in the next subsection we show that the bound-state wave functions 
of highly excited states agree with those of free fermions also for finite separations provided that the potential is much 
smaller than the energy, V <^ E. 

Denoting the 2-momentum of the current by P = (P'^ > 0, P^) the fermion loop amplitude in Fig. [t] is for a vector 
current 



The imaginary parts of the loop contribution for this as well as the scalar, pseudoscalar, and pseudovector currents 
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=11 



ImLsiP) 



pii pv , 



2m^ 



P2 ) ^p2(p2„4^2) 



^p2 _ 4j„2 



2VP2 - 4to2 



(vector : = ipl'^'ip), 
(scalar : js = V'V'), 
(pseudoscalar : — i^js'ip)^ 
(pseudovector : = ipj^j'^ip). 



(4.2) 
(4.3) 
(4.4) 
(4.5) 



The contribution of a bound state \n), with P^ |n) = P^J* \n), to the imaginary part of the vector current with 



P" > is 



ImL^J''(P) = Imz J d^xe'^P~P-^-^{0\j''{0)\n,t^O) {n,t^0\f{0)\0)9{x°) 
= 27TH'iP - PnmnO) \n) {n\fiO) |0) . 



The expressions (|3.1|) and (|3.7|) for an equal-mass bound state are 

dxidX2e*^"("i+"^)/2^(0, Xi)^ixi - X2)^0, X2) |0)^ , 



\n) 



^{x) = $o(a;) + $i(x)75 



2m ,Y 

—w 

a 



where 75 = 77. Using Tr (757^7") — 2e^'^ , where e" = 1, we get 

^(OIj'^(O) \n) = \n) = Tr [7^7"'i'(0)7"] 



4m 

'P2 



^Me^'^'P, 



(4.6) 

(4.7) 
(4.8) 

(4.9) 



Wc note that i?(0|j^(0) \n) P^ — Q as required by gauge invariance. The sum over inte rmediate states ^ 



includes an integral over the momentum P^ of the single bound state considered in (4.6). The (local) average over 
the energy P° of this bound-state contribution should be dual to the loop contribution (4.2), 



AM2 



(2^)2 



ImLJ^'^(P) ~ImL^'^(P). 



(4.10) 



Eq. (4.101 then gives 



Acc ording to (3.19) the bound-state separation AAf,^ — M^^i — Af,^ ~ 2tt at large masses [for states with $i(0) 7^ 0]. 

TT / A/2 

\<fiix = 0)f ~ - J (sector current). 



(4.11) 



Based on the form of the bound-state equation (3.9) we previously noted that $i(cr) is frame independent. Since 
a{x = 0) = is also frame independent this implies that $i(x — 0) cannot depend on P^. The duality relation 
(4.11) satisfies this constraint. 



We may determine the normalization analogously using currents with other Lorentz structures. Thus 

ii(0|Vi(0)V(0) \n) = Tr [7°*(0)7°] = 2$o(0) (scalar), (4.12) 

fl(0|Vi(0)75V'(0) |n) = Tr[757°$(0)7°] =-2$i(0) (pseudoscalar), (4.13) 



fl(0|Vi(0)757''V'(0) |n) = Tr[757^7°*(0)7°] =-^*i(0)P'' (pseudovector). (4.14) 



In D = 1 -I- 1 dimensions the pseudovector and vector currents are related, ip'y^'y^^ip = e'^''ip^u'4>, where 75 = 7''7"'^ and e^^ = 1. For 
completeness we anyhow give results for both. 
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The duality relation for the scalar current gives 



l*o(0)|' 




M2 - 4m2 
Ml 



(scalar current). 



(4.15) 



According to p.l7[ ) either $o(0) 
complements t 



or <i>i(0) = for any given bound state, so (4.15) is consistent with and 



re condition (4.11) for $i(0). The duality conditions on $i(0) obtained using the pseudoscalar and 



pseudovector currents agree with the vector current relation (4.11) 
For highly excited states we thus find that |$o(0)P 



asymptotic formula (3.16) we find that, up to an irrelevant phase. 



N = 



if $i(0) = [$o(0) = 0]. Using the 



(4.16) 



in (3.151. 



B. Wave function in momentum space at large M 



It is interesting to consider the duality indicated in Fig. [7| also more differentially: Does the wave function of highly 
excited bound states resemble the free fermion distribution given by the imaginary part of the loop, as expected in the 
parton model? For large masses, M ^ to, we may set m — 0. In the rest frame (P^ = 0) the free fermion momenta 
are k = iM(l,±l). 

The bound state (3.1 ) is in the rest frame 
dkidk2 



\M,Q) = j dxidx2 

The TO = spinors 
u{k) 

satisfy 



(27r)24|fci||fc2 

|fc|CT3 - kia2 ( 1 




u(A:i)e-*'=i^i5l^+«(fci)e''=i^idfe,J ^(x^-x^) [^/(A:2)e^'=^^^&fe,+«(fc2)e-''="^4J r ■ 

(4.17) 



v(k) 



\k\a-i -kia2 ( 
1 



(Tiu(fc), 



(4.18) 



u(k)u(k) = u{k)aiu{k) = v{k)v{k) ~ v{k)aiv(k) = 0, 

u{k)v{~k) = -v(k)u{-k) = -2k, 

u{k)aiv{—k) = —v(k)aiu{—k) = 2\k\. 

In terms of the Fourier transform of the wave function, 



ikx 



the state (4.17) becomes 
|A/,0) 



dk 



27r2|fc 



■{ [ - e(fc)a>o(fc) + $i(fc)]6l dlfc - [e{k)^o{k) + $i(fc)]d_fe 6^} |0)^ , 



(4.19) 



(4.20) 



(4.21) 



where e(fc) is the sign function. It appears to contain both positive {b\ c^lfc) and negative (d_fc 6^) energy modes. The 
wave function of the negative energy modes, however, turns out to vanish. 

is 



For TOi = TO2 — and — Q the bound-state equation (3.9) is 
idMx)^\{M~V)<^i{x), 
The solutions with a continuous derivative at x = are 



(4.22) 



$i(a;) + $o(a;) 



exp 



'-Mx+ -x^e{x) 



$i(a;) - $o(a;) = ?7w|exp ''-Mx ~ ''-x^e{x) 



(4.23) 
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with ^j{--x) = {—iy'^^ri^j{x) as in (3.171. The normahzation was determined (up to a phase) by the duahty 
conditions ( |4.11 1 and (4.15). In momentum space (4.20) we get 



$i(/c) + $o(fc) = / dxexp[-i{^M + k)x + ix^e{x)/8] ^ l{2n)^^^6{k + ^M), 



^i{k) - $o(fc) - ^v{2nf^^S{k ~ \M), 
where the approximation is vaHd in the large mass hmit, M ^ 1. Consequently 

e(fc)$o(fc) + $i(A:) = 0, 



and the bound state (4.21 ) reduces to 



|M,0) 



'2M 



7t 



7t 



M/2 "'-A//2 ' "-M/2 "'Af/2 



|0> 



R ' 



(4.24) 



(4.25) 



(4.26) 



with a momentum distribution of free fermions, in agreement with the parton model. The approximation made in 
(4.24) breaks down at large fermion separations where V{x) > M and the effects of confinement set in. Thus the 
fermions are approximately free only at shorter distances. 

An analogous study of duality may be made in a frame with non- vanishing CM momentum P. The parton state 
corresponding to (4.261 is then 



2M 



E/2+P/2 "-B/2+P/2 ' "-B/2+P/2 "'E/2+P/2 



|0), 



(4.27) 



V. ELECTROMAGNETIC FORM FACTORS 



A. Definition (mi 7^ 1x12) 



The electromagnetic current is 



f{z) = ^e/Vi/(^)7'^V'/(^) = e^^-j^(0)e 



(5.1) 



where e/ is the electric charge of flavor / and P = {P''\P^) is the generator of time and space translations. We 
consider the matrix element of j^(z) between bound states of the form (3.11, 



\A{Pa)) = / dxidx2 exp[liP^{xi+X2)]M0,xi)'i>A{xi-X2)M0,X2)\0)i^, 



(5.2) 



where we used the notation (3.33) for the 2x2 wave function ^a, whose structure is given by (3.5) and (3.7 1. Since 
the bound states are eigenstates of energy and momentum, P^ \ A{Pa)) = P^ \A{Pa)), the form factor can be expressed 
as 



Fj^^(z) = {B{P,)[f{z) \A{Pa)) = e'(^-^»)-(i?(n)|jnO) |A(P,)) 



(5.3) 



Due to the retarded boundary conditions, iIji{z) |0)^ = i/jliz) |0)^ = 0, only anticommutators between the fields of 
the current with those of the states contribute. In effect, the states \A) and {B\ are treated as the free \in) and 
{out\ stat es of standard perturbation theory. Here th e as ymptotic states are bound by the instantaneous Coulomb 
potential (1.10) arising from the boundary condition (1.6) on A" and have no O (e^) contributions. We expect that 
a perturbative expansion can be formulated analogously to the one for free asymptotic states. In the following we 
restrict ourselves to the lowest-order contribution. 



Since the potential is confining we consider only neutral bound states, and thus set ei = £2 = 1 in (5.1 ). Then the 
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current couples equally to both flavors 

2 



/=i ■' 



/2-j(yi+j/2)P6/2 



X fl(0|^|(0,y2)*L(yi-y2)7Vi(0,yi)^V^/(0,0)7^V'/(0,0)^i(0,a;i)^'A(a;i-a;2)V'2(0,X2)|0)j 



^.^-Pi)-/2 |Tr [vl/t^(x)7^7"vl/^(a;)] - ry^Ty^Tr [vl/B(a;)7%'^\I'^(a;)] }, (5.4) 

wh ere, in the second (/ = 2) term, we used (7'') 7° = 7^7^ (since 7° = 0-3, 7^ = 10-2) and \t'(— x) = ?75'*(x) according 
to K^. 



B. Gauge invariance {D = 3 + 1) 



Gauge invariance of the form factor (5.3) requires that 



G%{z)^dlF^/l''{z)=0 if = 1,2) 



(5.5) 



separately for the current of each flavor /. We shall show that (5.5) is a consequence of the bound-state equations 



satisfied by ^'^i and 5* b ■ Since the derivation is essentially independent of the masses and of the number of space-time 
dimensions we here consider mi ^ m2 and D = 3 -I- 1. Then 



Fji']^(z) = e*(^^-^")-^y"da;e'(^'-P")-/2Tr[*yx)7'y*^((r)], 



(5.6) 



Due to translation invariance we set 2; = in Gab{z) without loss of generality. The bound-state equations (1.91 for 
^Aix) and ^^a;) are 



iV, • {7*'7, ^a} - \Pa ■ [7%, + TOi7°*A - m2*A7° = (^a - V)-^A, 



(5.7) 



Multiplying the first equation by — ^'^ from the left, the second by from the right, and taking the trace of their 
sum gives 



-iV, -Tr (7%{*|j,*^}) +i(Pfc-P,)-Tr (7% [^Ij,*^]) - (i;^ - Sa)Tr (^^j^-^ 
Using [^I^^, = {vj/^, vj/^} - 2^1'^^'^ and multiplying both sides by exp[«(P(, - Pa) ■ x/2] we find 



(5.8) 



/2 Tr (7O7 [^\,-^a}) ] = e^(^-^»)-/2 Tr [l>l{x){fB ~ ^Ah'^'^Aix)] . (5.9) 



Integrating both sides over x the rhs. becomes — iG|^^(0) and the Ihs. vanishes (assuming that the integral over the 
oscillating wave functions is regularized as \x\ — 00, similarly as for plane waves). This proves the gauge condition 
(5.5) for f — 1. For / = 2 the gauge term corresponding to (5.6) is 



and the proof that it vanishes is analogous to the above. 



(5.10) 
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C. Form factor for mi = m2 {D = 1 + 1) 

The expression ( |5.4[ ) for the form factor simplifies in the case of equal masses, nii — m,2 — m. Since ip = the 
structure (3.7) of the wave function ^l* = $ becomes 



= <^>o{x) + $i(x)7°7^ + 2m$i(a;)— 7°7\ 

a 



(5.11) 



where n(a;) — {P'^ — V{x),P^) is the kinematical 2-momentum (3.4). The traces in (5.4) are now 



4m^ ~ „ 

1 + n, • Ub 



4to^ 



(5.12) 



where Ha = [P^ - V(x), -Pi) and Eqi = -1. 



The constraint (5.5) of gauge invariance implies that the form factor in Z? = 1 + 1 can be expressed as 



(5.13) 



where — —q^- Solving this for Fab{Q'^) with = 0, using Eq. (5.4) for the Z/is., and inserting the traces of (5.12), 
we obtain 



Fab{Q'') = -^ 



q^ Jq V 2 / L V aaCTb 



(5.14) 



The form factor vanishes unless "qarjh = —1- Therefore the factor in the square brackets could be taken to be 
antisymmetric in x, which allowed us to restrict the integration to positive x. 



VI. DIS AND PARTON DISTRIBUTIONS 



We consider the cross section for e(fci) + A[Pa) — >■ e(fc2) + B{Pi,) in the limit where Xbj = Q^/ (2Pa • q) is fixed, with 
q = ki — k2 and = —q^- The "inclusive" system is thus a discrete bound state B. The cross section is proportional 
to the square of the form factor Fab{Q^) defined in (5.13), with Mb oc Q. 



There are some peculiarities with DIS in D = 1 + 1 as compared to £> = 3 + 1: 

• The very concept of a "cross section" is related to transverse size. We may nevertheless define a Lorentz- invariant 
cross section by analogy to the usual case and then compare parton and bound state cross sections. 

• The virtual photon has a finite longitudinal ("loffe") coherence length in the target rest frame, L/ ~ Q^^ v/Q = 
l/(2maXBj)- In the absence of transverse dimensions DIS photons can be coherent on several partons at leading 
twist. The fractional momentum of a struck parton is kinematically constrained to be xbj- 

• In = 1 + 1 the scattering angle can be 6* = (forward) or 9 — tt (backward scattering). At the parton- 
level, forward (elastic) scattering implies = and thus is irrelevant for the Bjorken (Bj) limit. In backward 
scattering s ^ —t = , which is analogous to standard DIS in the Breit (or brick- wall) frame. 

• Since the coupling e has the dimension of energy in Z) = 1 + 1, the parton- level cross section will on dimensional 
grounds be suppressed (compared to D = 3 -I- 1) by a factor e^/Q*. 

• The backward scattering amplitude of elementary spin-^ fermions is proportional to the fcrmion masses. This 
further suppresses their scattering cross section. 
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A. The parton distribution 



The kinematics of DIS in 13 = 1 + 1 is discussed in Appendix |A] We find tliat the parton distribution may be 
expressed as 



XBj 



-\Q'Fab{Q' 



(6.1) 



where m is the mass of the target parton and the invariant form factor Fab{Q'^) is given in ( 5.14[ ) (for a neutral // 
state with mi — m2 = to). The mass of the inclusive system in the Bj limit is 



Ml^Q'^i—-! 



^XBj 



(6.2) 



In order to calculate the leading twist parton distribution for a neutral two-body state we analyze the expres- 
sion (5.14) at large and Mf,. We work in the Breit frame where q° = P° — P° = and = —Q is large. The basic 
expectation is that the Fourier phase in (5.14) limits the integration to x < l/Q, which is the Lorentz contracted 



equivalent of a finite (loffe) distance in the target rest frame. The integrand is O in this region, and the measure 
adds a factor of l/Q, so we obtain a contribution at leading order Fab{Q^) ^ Leading contributions can, 

however, arise also for larger (typically ~ Q") values of x, if the oscillations of the wave functions should cancel the 
Fourier phase such that a stationary phase arises. Such contributions are analyzed in detail in Appendix [B] and shown 
not to affect the leading result. 

It is convenient to introduce a rescaled variable 



(6.3) 



In the Bj limit, taking v = O (Q^) and using the expressions (A. 18) for the momenta, the variable a defined in (3.3) 
is of O (Q") for the target, 



2xBj 4 XBj 



while it is of O (Q^) for the final state. 



(6.4) 



a. = 0^(^-l)-^N + ix^-g^(^-l)-^.r.. 

\XBj / 2xBj 4 \XBj I XBj 



(6.5) 



Thus we may use the asymptotic forms (3.161 for ^B^Ph) 



-1 



The condition d^^^QB^y = 0) of (3.17) determines M\, such that cos[iTf, — rv? log(Tf,) + argr(l + irv?')\ = 1 
at w = (up to an irrelevant phase). The v dependence of the logarithm is of O (Q^^) and may be ignored. As 
the state is highly excited, we may use the normalization from (4.16). Thus 



*ob(t6) 



I 

2 



/ \ ^XBj ' 



\2xBj 

(—) 



(6.6) 



The expression (5.14) for the DIS form factor becomes 



Q^FABiVb = -) - -'iiV2^{l+Va) I dv smv\cos(-^)i^oA{Ta)~sm(-^)^iA{Ta)(l + ^^)]. (6.7) 

Jo L \2xBj/ \2xBjJ V XBjTaJi 
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For large v, Ta — > — oo and we may use^^ the asymptotic expressions (3.161 also for ^a- 

'2 N 



$0A(ra) ^ - - Ve2™= - 1 e-™' cos U{v) ~ 

{to) ~ W- — v/e2^™^^^^e-'^"%inf^(u) 

V TT TO V 



(y — ?> oo) 



(6.8) 



2xB-i 



where 



2 V 

^(v) = -TO log ^ ^ + argr(l 



and M = Mq is the target mass. The term in the square brackets in (6.7) behaves as 

2to2 



XBjTa 



2 N 

TT TO 



(6.9) 



cos/x(z;) (6.10) 



for t; — >■ cx). Because the integrand of (6.7) oscillates asymptotically, the integral does not converge. It can, 



however, be defined in the standard fashion by adding a "convergence factor" , e.g., a factor e in the integrand, 
and by taking e — > in the end. 

Recall that we assumed that Fab{Q^) only receives leading contributions from the region x ^ 1/Q in order 
to derive the result (6.7). The fact that the regularization procedure works suggests that this assumption was 



correct. We analyze this in detail in Appendix [B] 
• r]b = l 



The condition ^ob{v = 0) = of (3.17) determines Mf, such that sin[irfc — to^ log(Tf,) + argr(l + irn^)] = 1 at 
V = 0. Hence, 



sm , 
2 \2x 



\2xBj^ 



'2-K 



cos „ 
2 V2xi3,- 



f— )■ 



(6.11) 



The expression (5.14) of the DIS form factor becomes 



'FABiVb = +) =i -4jV27r(l - 77a) / dv sinv 



sm I :r-^^i$oA(Ta) + cos (— ^) $iA(ra) (l + 



2x 



2xb] 



XB]Ta 



(6.12) 



B. Numerical evaluation of the parton distribution 



Let us insert the explicit results (3.15) in the expression (6.7) and evaluate the integral numerically. We limit 
ourselves to the case i]a = —1% = +1- The case of opposite parities can be analyzed similarly. Since the integral does 
not converge we need to subtract the divergent term and treat it separately. We write 



^FABiVa = +)= h+h, 



(6.13) 



Recall that we already took Q — > oo, thus more precisely, the asymptotic expressions hold for Q 3> 2> 1. 



24 



Ii = -8i\/27r I dv smv cos (^^^)i$oA(Ta) - sin (^^^^$iA(rQ)(^l 

Ve^™' - le-™'cosAi(i;)l, (6.14) 



where 

2to2 

12 N ^ 
TT m 

/2 = Ve^™' - 1 e"''" / sinwcosAi(w). (6.15) 

™ Jo 

The divergence now only appears in the second integral, which may be calculated analytically. The standard regular- 
ization yields 

/2 = -16i^ v/e2-™' - 1 e"""" cosh (^^) cos + logXBj^ \T{l + im^)\ . (6.16) 

The first integral Ii can be calculated numerically^'^. The parton distribution can then be found by using the formula 

We show numerical results for the parton distribution in Fig. |8] The target wave function was chosen to be the 
one with lowest nonzero mass, which indeed has r]a = +1. We used Eq. (4.16) for the normalization of the target 
wave function, and chose two reference values m = 0.1 and m — 1.0 for the fermion masses. For m = 0.1 the target 
is highly relativistic, such that M ~ 1.78 ^ 2m, whereas for the second choice m = 1 we find that M ~ 2.70, i.e., 
a binding energy smaller than the constituent masses. This is reflected in the low-xsj behavior of f(xBj), which is 
qualitatively different in the two cases. The red curves show the XBj — ^ expansion of f{xBj), which is calculated in 
Appendix [Cj 

One can check numerically that the integral I2, which arises from asymptotic oscillations of the wave functions, 
dominates the result for small xbj and to. Therefore approximately 

XBjfixBj) - cos^ (mV2 + \ogXBj) (6.17) 



in this region. Inserting the result for the bound state mass at small to with n = 1 from (3.231 we find 

XBjf{xBj) ^ sin^ [to^ (log XBj + log TT - Ci(7r) + 7e)] ~ m'^ (log xbj + 1.648)^ . (6.18) 

Thus the contribution from the oscillations has a node at^^ relatively small xbj — 0.19 when to is small, and grows 
logarithmically ~ (logXBj)^ as xbj 0, until the distribution "saturates" at xbj ~ e^^^™^. 



VII. CONCLUDING REMARKS 



Analytical studies of bound state dynamics are usually based on summing Feynman diagrams, 6.17., using Dyson- 
Schwinger techniques |15l 116] . In the weak coupling limit the dynamics is non-relativistic and determined by the 
Schrodinger equation. Relativistic, confined states like hadrons may then emerge only when the coupling is strong. In 
D = 1 + 1 dimensions some all-orders, exact results have been obtained, notably for QED at zero fermion mass (the 
Schwinger model [TT]) and for QCD in the limit of a large number of colors, iVc —?> 00 (the 't Hooft model [H])- Their 
study have led to valuable insights ~ see, e.g., [T^VET] and references therein. Bosonisation of the massive Schwinger 
model furthermore allowed to obtain some approximate results in the strong coupling limit |22j . No analogous results 
have been found in D = 3 + 1 dimensions. Approximations based on a truncation of the Dyson-Schwinger equations 
have allowed analytical studies of hadrons in QCD [151 116] . which complement numerical results using lattice methods. 
Currently holographic approaches to QCD motivated by the AdS/CFT correspondence [23J[23] are under intense study 
as a means of obtaining results in the strong-coupling limit of the theory. 



This integral is still not absolutely convergent, which makes the numerical integration rather tricky. We found the best results by 
using the ExtrapolatingOscillatory method of NIntegrate in Mathematica. One can also introduce a cutoff at large v and use standard 
algorithms for the numerical integration. 

The approximation for the location of the node is poor since terms suppressed by XBj were neglected. 
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FIG. 8: The parton distribution of the target ground state. We plot the result for XBjf{xBj) as a function of xsj in linear 
scale (top) and in logarithmic scale for low XBj (bottom). We used m = 1 (left) and m = 0.1 (right) for the fermion masses. 
The dashed red curves show the asymptotic behavior of f(xBj) up to next-to-leading order as xsj — > [see Eq. (C.14|]. 



In this paper we explored a rather different approach to relativistic bound states in gauge field theory. It may be 
relevant provided the QCD coupling remains perturbative even in the long-distance regime governing hadron binding. 
Our approach is motivated by features of the data discussed in the Introduction, as well as by phenomenological and 
theoretical studies which indicate that as freezes at a moderate value in the infrared |25H30j . This possibility merits 
attention since it allows to bring the powerful techniques of perturbation theory to bear on bound state dynamics. 

In our scenario the confining potential arises from a boundary condition imposed on the solution of Gauss' law [5]. 
This gives an exactly linear potential even in D = 3 -I- 1 dimensions, with strength determined by a parameter 
A related to the boundary condition. Since A is independent of the potential is of leading order compared to the 
perturbative O (a^) interactions. Hence the bound state wave functions have no contributions from transverse gauge 
bosons at lowest order. 

Quark-antiquark pairs will be created even in the lowest order potential, leading to hadron Fock states with 
any number of fermions. This is also true for the Dirac states of an electron in an external potential. The Dirac 
wave function, which has a single electron degree of freedom, is obtained using retarded boundary conditions for 
the electron propagation. Hence the norm of the Dirac wave function gives the inclusive electron density. In our 
framework mesons are analogously (at lowest order in Us) described by ostensibly two- particle {qq) wave functions, 
which implicitly include such quark-antiquark pair production effects. The use of retarded boundary conditions for 
the quarks is allowed in the external legs of S'-matrix elements since they are amputated and put on-shell. In analogy 
to the standard perturbative expansion of scattering amplitudes based on the free in and out states of elementary 
fields we envisage an expansion based on the Born-level bound states formed by the linear potential. 

Previously [7] we verified that the // Born states are Poincare covariant in D = 1 -I- 1, as expected because the 
linear potential arises from a boundary condition which is compatible with the field equations of motion. The wave 
functions turned out to depend on the separation x between the fermions through the frame-invariant square IV^ — a 
of the "kinematical momentum" li = {E —V,P), where {E,P) is the 2-momentum of the bound state and V{x) is 
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the linear potential. An earlier study [5] indicated that Poincare invariance is preserved also in D = 3 + 1. 

In this paper we found the analytic expressions of the Born level // wave functions in Z) = 1 + 1 and studied their 
properties. Their norm tends to a constant at large x. The nearly non-relativistic case shown in Fig. |4] makes it 
clear that the constant norm reflects fermion pair production at large values of the potential. The norm of the Dirac 
wave function behaves similarly (Fig. |2|, which in that case implies a continuous Dirac mass spectrum^^ [IO] - [l3] . The 
// wave function is, however, generally singular at the value of x where the "kinematical mass" vanishes, iP = 0. 
Requiring the wave function to be regular gives a discrete spectrum. 

We found that the wave functions of highly excited bound states agree with parton model expectations in the range 
of X where the potential is small compared to the bound state mass: Only fermions and antifermions of positive energy 
contribute to the bound state, having the momenta of nearly free particles. Duality allows to determine the overall 
normalization of the wave functions through the condition that the (average) contribution of the bound states to the 
imaginary part of current propagators agrees with that of free fermions. The duality relation works in all frames and 
for all currents. 

A manifestation of the virtual fermion pairs in the bound states is provided by the parton distributions measured 
by deep inelastic lepton scattering. For relativistic states the parton distribution grows as XBj — > 0, qualitatively in 
agreement with data on sea quarks (although the present calculation is in one space dimension only). The use of 
retarded boundary conditions allows a simple "valence" description of bound states, while preserving its multiparticle 
features. 
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Appendix A: DIS in D = 1 -|- 1 
1. Units and kinematics in Z) = 1 + 1 

The dimension of the Lagrangian is [C] = E"^ in energy units. Hence the fermion and photon fields have [ijj] — E^^^ 
and [A^] = -E", while the electron charge^^ [e] — E^. With standard spinor normalizations uu = —vv = 2m the 
fermion operators have [b] = [d] = E^ and satisfy {b{k^),b\p^)} — 2Ep2-K5{k^ — P^)- The states have the same 
dimensions as their creation operators, [|e~)] = E^ . For our bound-state definition 
Consequently the invariant form factor defined in ( |5.13| is dimensionless, [Fab{Q'^)] = E". 

The 2—^2 scattering amplitudes 

A{e^j. ^ e^i) = {e^i\T \efi) ^ M{efi ^ e^i){2TTfS'^{P, - Pf) (A.l) 
have [A] = E^ and [A^] = E^ . Defining the 2-^2 "cross section" in analogy to D = 3 + 1 as 

a3cat(eM->e/i) = ^ j (2^^^^l^l'(2'^)''^'(^'i +^^2 - (A.2) 



5^ this implies [$] = E'^. 



The Dirac spectrum is continuous for almost all potentials in both two and four dimensions. Curiously, textbooks rarely mention this 
fact, even though the exceptional case of the 1/r potential in D = 3 + 1 is treated in detail. 
^® For clarity we display the charge e explicitly in this Appendix, taking V{x) = ^e^\x\. 
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gives [cTscat] = as expected. The phase space factor is 

dpi 



1 \M\^ 
2s 4e^,p5'p^ 



where e*^^ = 1, and we used 



d{E^+E^) _p\ p\-Pl _ 1 10 „ 1 _ e^.p^p- 

.1 ^ + ^ p. p. J - 



E1E2 



dp\ El E2 E1E2 

The invariant kinematic factor may be evaluated in the Breit frame where, for Q ^ mi, m2, 

Pl = f(l,l), P2-|(l,-l), 



givmg 



4Q4 



for ^ 00 with fixed mi, m2. 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



We may evaluate the electron vertex as 



m „ + m I 



^{Ei +m){E2 + m) 



{El + m)k!^ + {E2 + m,)fcf + i.g^°(fci - ^2)" = e'^'^g,y(Q), 



(A.7) 



where q = ki — k2. For backward scattering, taking the incoming electron momentum kl < and thus ^2 > 
may evaluate explicitly as 

1 /7 T /7 77 T (-El + m) (£'2 - m) - f^;! - m) (E'a + m,) q" , . . 

= ^[Ei+m){E2-m)-^{Ei-m){E2+m)= ^ \ ^ ^i^^^ ; ^ = -2m^, (A.8) 

^(£■1 + "i)(£;2 - m) + ^(£;i - m)(£;2 + m) Q 



which implies 



V{Q) = 



2m 



for the invariant part of the vertex in (A.7). Using 



the backward scattering amplitude for e/i — > e/i is 



Using this in the expression (A.6) for the cross section at large Q gives 



CTscat (e/i e/i) 4m^m^ — , 



(A.9) 



(A.IO) 



(A.ll) 



(A.12) 



which is suppressed by m^mj^/Q compared to the "scaling" behavior ^ e /Q . 

We may compare the fermion cross section with that for elementary scalars, which do not have a mass suppression. 
The scalar vertex is 



y/ = (fci + fc2)^ = e'^"g.v;(g). 



(A.13) 
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We can determine Vs{Q) from the space component (taking kl < 0) 



implying 



Consequently the scalar 2 — >• 2 amphtude analogous to (A. 11) is of O (Q"), giving a scaling cross section in (A.6|, 
CTscat oc e^'/Q^. 

2. e(fci)A(Pa) e{k2)B{Pb) in the Bj limit 

The Bj limit is defined as usual, 

Q^ = -iPb-Paf^^ with XBj = T^^ fixed. (A.16) 
The mass Ma of the target A is kept fixed, while the mass of the produced bound state grows with Q, 

m2 = (P, + g)2 ^ g2 _ iV (A.17) 

In the Breit frame, 

g=(0,-Q), fci = iQ(l,-l), fc2 = ^Q(l,l), 

= ^(1,1), Pb = -^(l,l-2a;B,). (A.18) 

2x3 j 2xbj 



Using the expression (5.13) for the bound state vertex we find 



M{eA^eB) = -e^^FAB{Q'')- (A.19) 



The kinematic factors in the expression (A. 3) of the cross section are 



s 



Q 



2 



Hence the DIS cross section is 



and e^^ki^P^ = \Q\ (A.20) 



a.cat(eA ^ eB) = e'^^XBAFABiQ^)? ^ Sttg^^, (A.21) 



b 



where the second equality implies an average over the bound state peaks, whose separation in Af^ is 27re^ according 
to (3.19). Converting dM"^ = —Q'^dxBj / XBj^ we find the parton distribution f{xBj) of the target as 



= ^1 7^\Pab{Q^)? = ^sci,t{ey,^ en)f[xBj), 

dXBj iTTXBjQ 



1 1 

Sne^m? xbj 

where m is the mass of the struck fermion in the target A. 



/(^s,) = ^^:^ — \Q'FABiQ^)\', (A.22) 
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Appendix B: Details on the Q'^ — > oo limit 



The leading O (Q^^) result ( [6?7t, (|6.12[ ) for the DIS form factor Fab{Q'^) in the Bj hmit was found by calculating 
the contribution to the integral ( 5.14^ for x ^ l/Q explicitly. Leading contributions can, in principle, arise also for 
larger values of a;, if the oscillations of the wave functions cancel the Fourier phase such that a stationary phase arises. 
In this Appendix we check that such extra contributions are absent. 



The form factor was given in (5.141 and becomes 



. 1 - rjarib 



Q 



Qx 



4to^ 



n, • n, 



(B.l) 



in the Breit frame. It is necessary to check the behavior of the expression in the square brackets of ( |B.l I by using the 
asymptotic expressions for $j. The wave functions depend on x through the variables cJa.b- The variable corresponding 
to the final state is 



\XBj ' ^XBj 4 



(B.2) 



which is large for a; > except very close to the roots 

x±= Q 







xbj ^ 


KXBj J_ 



(B.3) 



The asymptotic expansion of $3 can thus be used unless |x — a;±| < l/Q- The target wave function depends on 

Q , , 1 



aa^M^--^\x\ + -x\ 
2xB, 4 



(B.4) 



The asymptotic formulas for $^ can be used unless x < 1/Q (which was already discussed in the main text) or 
\x~2Q/xBj \ < l/Q. 

We shall discuss the asymptotics of the wave functions as a — > — cx). To next-to- leading order we find [compare 
to (plel)] 



$i((t) 



277771^ 



1 e^''"' <^ sin [f - TO^ log(-(T) + argr(l + im^)] 



$0(0-) 



V TT m 

TO** COS [| — 777? log(— cr) + arg r(l + zto^)] + sin — log(— cr) + arg r(l + im^)] 

a 

-i\ - — Ve2™' - le""™' \ cos [f - log(-cr) + arg Tfl + im^)] 
V TT m 

sin [| — log(— cr) + argr(l + im^)] + cos [| — log(— cr) + argr(l + im^)] 



(B.5) 



+ 0(a- 



Also, the expression appearing in the last term in the square brackets of ( |B.1[ ) can be written as 

il^fr TT -9 2 g^a + Q'b + (1 - XBjfQ"^ IxBj'^ 



(B.6) 



Notice that the apparent singularities of this term as aa,b are regularized by the zeroes of in (B.l). 

We start by discussing the contributions from the regions where x > O {Q'^) and the asymptotic formulas (B.5) 
work. Then the aa,b are O (Q) or larger. Thus the factor (B.6) is suppressed by l/Q, at least. Neglecting this factor, 
let us first consider the leading terms of the asymptotic expansion (B.5). The leading terms in (B.l) then combine, 
through sin(Q!a) sin(a6) + cos(Q;a) cos(Q!fc) — cos(Q;a — af,), to give 



FABiQ^) 



. 1 - VaVb 



-I -i — '— I dx sin f — — ^ cos 

Q .1 \ 2 J 



— 7;—+m log 



(Jaix) 



(B.7) 
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where, for Pi 







PS 



E, 



(B., 



is fixed and x independent. The remaining x dependence is logarithmic and cannot cancel the rapidly oscillating 
Fourier phase^^. Hence no stationary phase can arise from the leading behavior of the first two terms in the integrand, 
which suggests that the integral is limited to xQ of O (1). Notice also that we used the approximation of (6.4| and (6.5) 
in the main text, which is valid when x The variation of (B.7) due to this approximation is 



1 
Q 



dx 



'Qx\ X 



1 



(B.9) 



and thus indeed subleading. 

Let us then discuss the inclusion of the next-to- leading terms in (B.5 1 or the term proportional to the factor ( |B.6| ). 
These terms are suppressed by at least one power of aa or (Tb, i.e., a power of l/Q. As the integral in (B.l) is only 
multiplied by l/Q, they could still contribute at O {Q~^) to the form factor if a suitable stationary phase arises. A 
stationary phase at a generic x ^ Q is not relevant, because then (Tq ~ Q^, which leads to a suppression of at least 
Q~^. As we shall see below, the stationary phases only occur for a range of x with length O (Q^) , so additional powers 
of Q cannot arise from the integration. A stationary phase at a; ~ would, however, lead to an extra contribution 
to the form factor at leading order^^. 

Using the expressions (B.5I and (B.6I to expand the factor in the square brackets of (B.l) to next-to-leading order 
in 1/tT, we observe that all possible combinations of the Fourier phase and the phases of ^a.b appear. Neglecting 
constant and slowly varying factors, they are proportional to Qx ± aa i <Jb^ where the signs of and can be 
different. If they are, however, the situation is as in the case of the leading order analysis in (B.7), and no stationary 
phases arise. Therefore, we discuss only the phases Qxi [a a + which were absent at leading order. The locations 
of the stationary phases are found by solving 

[Qx ± K + a^)] = 0. (B.IO) 



dx 



The solutions are given by 



Q 



± 1 



XBj 



x±, 



(B.ll) 



i.e., they occur at generic O (Q) values. Near these points the phases behave as ^ \{x — x±)'^, such that the stationary 
phases are limited to regions having lengths ~ Q^, as expected. We conclude that the next-to- leading terms do not 
contribute to the form factor at leading order^^. 

Finally, let us discuss the contributions to Fab{Q^) from the regions where the asymptotic formulas do not hold 
for either of the wave functions '^a.b [so that da = O (Q") ot cti, — O The case x ~ l/Q was discussed in the 

main text. The other regions are the neighborhoods of x = x± in ( B.3[ ) and x — 2Q/xBj, respectively. Let us discuss, 
for definiteness, x = 2QJxBj. Near this point the asymptotic formulas for $b are valid, and we can write down an 
expression similar to (6.7| and (6.12), 



FabIQ^) 



1 

Q 



Qx 



dx sin 

l\x-2Q/xBj\r^l/Q 

where • • • stands for a function of xQ which is nontrivial but regular within the region of integration, 
integration variable we obtain 



Fab{Q'') 



1 

Q 



dx sin 



kl~i/Q 



^QX_ r 

\ 2 XBj ) ^ 



1 

Q^ 



dv sin 



■ ( ^^^\\ 

ni u H 

V XBj^ 



(B.12) 

Shifting the 

(B.13) 



The contribution from this region has thus the leading power behavior ~ 1/Q^, but also involves the large phase 
factor Q"^ /xBj. We interpret that the basically arbitrary phase averages the result to zero. Analogous results are 
found in the neighborhoods oi x = x±. 



Notice that this argument does not work when cFa or (T;, is close to zero, and the logarithmic term varies rapidly. This happens, however, 
only in the regions where the asymptotic expansions are not reliable to start with. 
Stationary phases near the roots x± of ( |B.3| or near x = 2Q/xBj would also be special. 

Notice that if 1 — XBj ~ so t hat t he final state has a finite mass, one of the stationary phases moves to x ^ Q*^ signaling the 

breakdown of the results \(y.7\ and ( |6.12| . 
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Appendix C: Asymptotics of the parton distribution at small xbj 



It is possible to calculate the xsj limit of the parton distribution analytically. We again assume that rja 
-rji, = +1 and start by separating the two contributions in Eq. (6.71 as 

Q^FAB{Va = +) = J1 + J2, 



Ji = — 8iv27r / dv sinv 
Jo 

16^m2^/2^ f°° , 

J2 ~ / dv sin 

XBj Jo 



cos 



/ V 



m V sm I , 

\2xBj'' Ta 



(C.l) 
(C.2) 
(C.3) 



where Ta in (6.4| is a function of v/xsj- It would seem that the expansions can be calculated by substituting v — )• XBjV 
in each integral and then developing the factor sm{xBjv) as a series at xbj — 0. This, however, leads to integrals that 
are divergent for v — )■ 00. Instead we expand the wave functions up to next-to-leading order for Ta — )■ —00, 



^l{ra) 



W2tt 
2 AS 

Phi 



2 



>os - ^) + {m^ - M^) sin {.(v) - 

(^(^)-2^) 



XBl 
V 



XBj 
V 



o 



XBj 



o 



XBj' 



(C.4) 



(C.5) 



where is given in (6.91, and we used N from (4.161 for the target wave function 



Let us discuss the integral Ji first. Using the expansions of the wave functions wc find 
J;^ = -8i7re^™^ y sin t;| cos -t- cos + ^m^ - sin^(w) — ^ 



( XBj"' 



\ f2 



(C.6) 



The integral arising from the first term in the wavy brackets was already evaluated in Eq. (6.16 1. The second term can 



also be calculated analytically. The first term in the square brackets contains a rapidly oscillating phase as XBj — 0. 
Hence its leading contribution arises from the region v ^ xbj, and is of O (^{xsj)^)- The dominant contributions to 
the second term have v ^ (xbj)^, and therefore the result is of O (xBj)- The O (^xbj^ /v^) term converges fast enough 
both as u — ?> and as v 00 for us to develop the sine factor as a series at v = and see that this contribution is of 
O ((a;_Bj)^). Altogether we get 



Ji = —Sine 



\T{l + im^)\ 



cosh 



cos 



' log XBj 



+XBj{m'^ - M^] sinh 



A I 

sm I + logXBj 



0{{xBj) 



(C.7) 



The integral J2 can be analyzed similarly. There are, however, complications due to the explicit factor of 1/xBj 
appearing in the coefficient in (C.3). First, we need to study the terms of the asymptotic expansion of the integrand 
up to O (xBj^/w^)- This is necessary in order to determine all O {xbj) contributions to the form factor from the 
region of asymptotically high v ^ (xbj)'^- Second, also the region with small v ^ XBj, where the wave function $1^ 
cannot be estimated in terms of elementary functions, contributes to the form factor at O (xBj)- It is hard to find a 
closed form expression for this contribution, but we will write it down as an integral below. 

Let us start with the contributions arising from the region with v ~ (xbj)'^- Developing the integrand at d — >■ 00 
gives 



J2 = Sim^TTe ^ I dv smv< 



-) 

XBj J 



cos ) 



(C.8) 



+(m^ — M^) cos i^Jiiv) ^ I — cos ^liv) — vr? sin ( /-ifi;) —\ + ™^ sin [liv] 

L V XBj' ^ a;Ri/ 



v_\ 

XBj' 
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The integral over the O (1/w) and O {xsj/v^') terms can be done analytically^" . This contribution is given by 



J'. 



2A 



8ie 7re 



|r(l + im ) I <^ sinh 



I sinh (^) cos + rn' hgxs,) ~ xs.e-^^'/' [{2m' - Al') sin (^) 



+w?(rr? - M2)cos + XBj{m' ~ M^) cosh (^^^^ sin +m'logXBj^ \ +0{{xbj)') ■ (C.9) 

The remaining contribution comes from small v ^ xsj It can be isolated by subtracting from the integrand its 
leading terms given in (C.8), which leads to 



SiTTin 



2 roc 



dv sinu-^ 



XBj Jo 



'iXBj I Ta 



-) 

XBjy 



COS ii{v) 



XBi 
V 



(C.IO) 



COS 



( ] — COS /iff) — m' sin ( /ifw) | + sin fj,{v) 

\ XBy V a;R,/ 



After scaling the integration variable by xbj the integral reads 

^ . ^V\^iA{Ta = AP -V) 



J'. 



2B 



SiTTm 



dv sm{xBjv)<^ 2 



cos 



X 



Bj 



COS (/i(w) — w) — cos ) — sin (/i(w) — t') + sin /!(?;) 



— u) — cos /2(u) 
1 1 



(C.ll) 



where il{v) — ^{xBjv) — M^/2 — m' logu, and therefore the whole expression in the curly brackets is independent of 
XBj- The leading term as xbj — >■ is then found by expanding the factor sm{xBjv) at small xbj and v, 

J2B = iC{m, M)xBj + O {{xBj)^) , 

where 



(C.12) 



C(m,M) = 8W 



dv< 2 



2 . fv\ 

—V sm - 

TT V2/ 



NP - V 



; (/2(f) — w) — cos/i(u) 



-(m^ — M^)e cos (/i(w) — i;) — cos/i(w) — sin (/i(w) — + sin/i(i') 



(C.13) 



Collecting the results, 
Q'^FAB{Va = +) = Ji + J2A + J: 



2B 



8i7re~^™^/^ |r(l + im')\ | cos + logxsj) - XBj{m'' - sin + m' \ogXBj 

(2m2-M2)sin(^ — j + 777,^(7712 _ 7v/2) J I + ^(7(771, M )2;sj + O ((xbj)^) . (C.14) 



The XBj^ expansion for XBjf{xBj) obtained by using this formula in (6.11 is shown by the dashed red curves 
in Fig. [sj Notice that the contribution from the asymptotic oscillations is suppressed by the factor e"'^'^™^/^ \^ ^j^g 
non-relativistic limit 777 — )■ 00. Such a suppression factor is absent in C{m,M) of ( C.13| , which mostly arises from 
small V ^ XBj, j.e, the region where the wave function $1^ is non- vanishing in the nonrelativistic limit. 



[1] J. Beringer et al. [Particle Data Group Collaboration], Phys. Rev. D 86 (2012) 010001. 



Notice that the integral over the O {xbj/v'^) term is convergent at 1; — > despite the factor of thanks to cancellations in the 

numerator. 

To be precise, Eq. | |C.9| already contains some O (xbj) contributions from the region v ^ XBj- These are the terms without logarithmic 
phases. 
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